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Fermi-Bose Correspondence at Finite Temperature 

Girish S. Setlur 
Department of Physics and Materials Research Laboratory, 
University of Illinois at Urbana- Champaign , Urbana II 61801 
The correspondence between fermi-sea/bose-condensate displacements and the number- 
conserving product of two fermi/bose fields is generalised to finite temperatures. It is shown 
that the straightforward generalisation that involves making the sea-bosons participate in 



the thermodynamic averaging does not work out. A remedy is found which involves treating 

OS 

the sea-displacement bosons at zero temperature while all finite temperature effects have 
to be lumped into the coefficients. We also show that all the finite temperature dynamical 

P 

four-point and six-point functions come out correctly as do some of the commutation rules 
involving the fermion/boson product that go through unscathed. It is also shown that this 
unusual prescription of leaving out the sea-bosons from the thermodynamic averaging does 
^ ■ in fact reproduce the correct RPA dielectric function at finite temperature. This article is 

OO ' not however, critic-proof. 



I. INTRODUCTION 



The attempts of recent preprints have been to construct a nonperturbative theory of interacting 

S" 



fermions/bosons. In order for this to have practical consequences, we would like to be able to generalise these 
ideas to finite temperatures where one might be able to study important physical phenomena such as phase 



transitions. The generalisation to finite temperatures is not as straightforward as it seems. The obvious gen- 
eralisation that involves making the sea-displacement bosons participate in the thermodynamic averaging 

•i—t ...... 

while retaining the original formula that relates the fermion/boson product to the sea-displacements does 

not work out. This is explained in detail and a remedy is found. The remedy involves making the following 

unusual prescription. The sea-displacement bosons are assumed to remain at zero temperature and do not 

participate in the thermodynamic averaging. Rather, the finite temperature effects have to be lumped into 

the coefficients thereby altering the relation that connects the fermion/boson product and sea-displacements. 

It must be made sure that in the bargain, commutaion rules involving the number-conserving products are 

not altered. All this is shown explicitly in the next few sections. 



II. FERMION PRODUCT AND SEA-DISPLACEMENT CORRESPONDENCE 

To motivate the development of the finite temperature correspondence let us first try to argue that the 
naive generalisation that seems very natural does not work out. One would expect from our earlier preprintB 
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that at finite temperature, the correlation functions in the bose language are merely thermodynamic ex- 
pectation values of the corresponding dynamical quantities. This does not seem to work out as this simple 
demonstration shows. Take for example, the momentum distribution at finite temerature, 

(4c k ) = n F (k) +^A k _ qi/2 (-qi)(aj c _ qi/2 (qi)a k _ qi/2 (qi)) 
qi 

- A k+ qi /2(-qi)<a k+qi/2 (qiK+ qi /2(qi)) (i) 

qi 

Since, 

Assume that the n F (k) are all evaluated at zero temperature as in the previous preprints. Taken at face 
value, one is obliged to compute the thermodynamic expectation values of the bose occupation probabilities 
assuming the chemical potential for the bosons is zero as the total number of fermions commutes with the 
sea-displacements (this means that we are allowed to create and destroy any number of bosons). Now such 
a calculation yields an infinite answer for (c^c^) as the sum over all qi diverges (is proportional to the 
volume of the system) . This means that we are no longer allowed to treat the sea-displacements bosons as 
participating in the thermodynamic averaging. Rather, they merely provide the correct time-evolution of 
the product of fermi fields, the thermodynamic aspect has to be absorbed into the coefficients n F (k) and 
A k (q), as we shall soon see. For this, let us repostulate a form for the product of fermi fields and argue 
why this should be as it is. The only guide is that all the finite temperature dynamical correlation functions 
involving the number conserving object c k +q/2 Ck -q/ 2 should be correctly reproduced. The commutation 
rules involving these number conserving products should not be damaged in the bargain. Let us start with 
the simplest case namely, ( c k + q /2 Ck -q/ 2 ) • We know what the answer should be, that is, 

( C I+q/2 C k-q/2) = ^q=0^F,/3(k) (3) 

where, 

n Fi/3 (k) = [e:rp(/3(e k - + l]" 1 (4) 
For this to happen we have to do the following, 

(1) Treat the bosons as before assuming that they are always at zero temperature and with zero chemical 
potential. 

(2) Assume that all finite temperature effects are lumped into the coefficients. 

(3) Fix the coefficients as before by demanding that the finite temperature dynamical moments of the 
number-conserving products of fermi fields come out right. In order to cut a long story short, let us merely 
postulate the form of the product and later on argue why this should be as it is. The following formulas are 
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true ONLY for temperatures above zero. At temperatures exactly equal to zero the answers have already 
been given elsewhereS. For q^Owe have, 



=!+ q /2 c k-q/2 = (W^)[A£(q)a k (-q)+A£(-q)4(q) 



+Ti(k,q) a k+ q /2- qi /2(qi) a k-qi/2(qi -q) 

- r 2 (k, q) ^ a L q / 2 +q 1 /2( c ii) a k+qi/2(qi - q) (5) 

qi#o 



Here, 



Ti(k,q) = J (1 - n Fji8 (k + q/2) V (1 - nir ,3 (k- q/2)) 



T 2 (k, q) = Jn Ftf3 Qi + q/2)n Fj/3 (k - q/2) 



A k(q) = V n ^( k + - nF ^ k - ^/ 2 )) ( 6 ) 

It can be shown quite easily that all the dynamical four and six-point functions of the object c k_|_q/2 c k-q/2 
are recovered exactly (See Appendix). For q = we have to find another expression in order for the kinetic 
energy operator to come out right. For this let us postulate a form and then compute the coefficients. 

4 c k = UfAK + Y S , i(k,qi)ak_ qi/2 (qi)a k -q 1 /2(qi) 
qi#o 

- Y 5 '2(k,qi)aj c+qi/2 (q 1 )a k+qi/2 (qi) (7) 
qi#o 

The coefficients Si and S2 here have to be chosen so as to recover the following form of the kinetic energy 
operator. The latter is necessary to ensure that the time-evolution of the various operators come out right. 
The correct form of the kinetic energy operator is therefore (taking a cue from the zero temperature caseEl), 

x = E(^)4(q)«k(q) (8) 

k,q 

This means, 

ek+ q /2^i( k + q/2, q) - £k- q /2S 2 (k - q/2, q) = (-5) (9) 
The obvious solution to this is, 

Sl(k + q/2,q) = S 2 (k-q/2,q) = l (10) 
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Thus we may rewrite the above formula for the occupation number operator as, 

c^Ck = n FiP (k) + aL qi/2 ( c li) a k-qi/2(qi) 

~ E a I+q 1 /2( C H) a k+qi/2(qi) (11) 

qi#o 

It is straightforward to show that these formulas reproduce the correct commutaion rules described below 
(q^O and q ^ 0) at least in so far as the terms linear in the sea-displacements are concerened. The claim is 
that practical results such as the answer for the dynamical current-current and density-density correlations 
should have the right form even though many of the nuances may not be right. Think about it, we have 
a whole class of exactly solvable models in any number of spatial dimensions and are able to recover RPA 
and all the other well-known results in terms of quadratures, a feature that is likely to persist even after 
introducing more complicated physical processess such as interaction with phonons. 

i C l' C k' ' C !+q/2 C k-q/2] = <V C k -q/2<5 k ' k+ q/2 - C k +q/2 C k' V ,k-q/2 (l 2 ) 

and, 

[ C lc+q/2 C k-q/2, C [' +q < /2 C k' -q' / 2 ] = 4+q/2 C k' -q'/2 <*k-q/2,k' +q'/2 ~ C k' +q' / 2 C k-q/2 (5 k' -q'/2,k+q/2 ( X3 ) 



A. Finite Temperature RPA 

In this subsection, we demonstrate that the RPA dielectric function is recovered exactly by selectively 
retaining parts of the coulomb interaction that lead to RPA. We know that the kinetic energy in the bose 
language is given by, 

^™ = E(^M(qH(q) (14) 

k,q 



For this let us choose, 



where, 



H * = E (15) 

q#0 



Pq = Et A k(q) fl k(-q) + A£(-q)a£(q)] (16) 

k 

From this it may be shown that the RPA dielectric function is recovered as the following demonstration 
shows. Assume that a weak time- varying external perturbation is applied as shown below, 

H ext = ^(f/ext(q,i) + f/; xt (-q,t))p q (17) 

q#0 
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here. 

U ext (r,t) = U e**-™* (18) 



Let us now write down the equations of motion for the varios bose fields 
,|<4(q)>=^ k (q)(4(q)) + (^ 



i -(aL(q)} = "k(q)<4(q)> + )A£(-q) (-q)(a^ (q)) + Af, (q)«t (-q))] 



k' 

+ (^ex t (q,*) + ^ t (-q,*))Af(-q) (19) 
-i^<a£(-q)) = ^(-qXa^-q)) + (^)Af(q) £[A£ (-q)<a£, (q)) + Af, (q)<a£ (-q))] 

k' 

+ (C^t(q,*) + t/ e *, t (-q,t)K(q) (20) 
Now let us decompose the expectation values as follows, 

(4(q)> = E/ ert (q,t)Ck(q) + ^(-q,t)I> k (q) (21) 

(4 f (-q)) = U: xt (- q ,t)Cl(-q) + U ext (d,t)Dt(-d) (22) 

(a k f (-q)) = U: xt (- q ,t)Cl(-q) + U ext (^t)Di(- q ) (23) 

The coefficients C k (q) and q) satisfy, 

w C k (q) - c k (q)C k (q) + (^K(-q) £[A£ (-q)C k , (q) + A* (q)D*, (-q)] + A^(-q) (24) 

k' 

-^(-q) = ^k(-qM(-q) + (^)A^(q) £X , (q)D*, (-q) + A£ , (-q)C k , (q)] + A£(q) (25) 

k ' 

Now, effective potential may be written as, 

Ue//(q,t) = ^ext(q,*) + (y)<P-q>'^ext(q,t) (26) 

here, 

(p- q ) - f/ext(q, *)<p- q )' + U: xt (-q, i)(p_ q }" (27) 

Using the fact tthat, 

(P-q)' = E A k(-<l) C k(q) + E A k (q)^(-q) (28) 

k k 

Solving these equations and using the fact that the dielectric function is just the ratio of the external divided 
by the effective potential we get, 



Vq n F ^(k + q/2) - n F:(3 (k - q/2) 
V ^ 

k 

Which is, lo and behold, the RPA of Bohm and Pines 



(q, tt) = 1 + f E M 47 _Z = 47 ^ (29) 

k U m 



III. BOSE PRODUCT AND CONDENSATE-DISPLACEMENT CORRESPONDENCE 



The generalisation to bose systems is entirely analogous. The answers will merely be quoted here, in fact 
it is easier to verify the various assertions in the bose case than it is in the fermi case. Again these formulas 
are valid strictly for temperatures above zero. At exctly zero temperature, the answers have been given 
elsewherei. 



&k&k = n B ,/3(k) + d k-q 1 / 2 ( ( li) d k-q 1 /2(q^ 



-qi/2 

qi#o 



2 4+q 1 /2( C H) d k+q 1 /2(qi) (30) 

qi#o 



and for q^O, 



& U/2^-q/2 = (^)[A^(q)d k (-q) + A^(-q)4(q)] 

+T 1 (k,q) J2 d lk+q/2-q 1 /2( f ll) d k-qi/2(qi - <j) 

qi#o 

+ ^2(k,q) ^ 4-q/2+q 1 /2(qi) d k+q 1 /2(qi ~ q) (31) 

qi#o 



and the coefficients are similarly given as, 



T!(k,q) = J(l + n s , /3 (k + q/2)),/(l + n s , /3 (k-q/2)) (32) 



T 2 (k, q) = Jn B ,p(k + q/2Wn Bi/3 (k - q/2) (33) 



and, 



and, 



A^' B (q) = V'4( k + <l/2)(l + 4(k - q/2)) (34) 



n£(k) = : (35) 

There are two points worth noting here, one is the change of sign in (1 + rig(k — q/2)) in the fermi case it 
was (1 — rip(k — q/2)), the reason should be obvious. The other is the overall change in sign in T2, this is 
necessary to ensure that the correct six point function is recovered. The proof that all the six-point functions 
come out right may be seen in Appendix. 
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IV. APPENDIX 



Assume that in the fermi case the part quadratic in the sea-displacements is given by, 



'k+q/2 C k-q/2|Quarf 



£ r£ ; £(k,q)a ki ( qi K 2 (q 2 ) 



ki,k 2 ,qi,q2 



From this we may write, 



1 - (4+q/2 Ck -q/ 2C l'+q72 C k'-q'/2 C l"+ q "/2 C k"-q"/2) 



[(1 - n FJj (k - q/2))(l - n F ^(k - q /2))n F ,p{1i + q/2)5 k+q/2;k »_ q » /2 5 k _ q/2;k ' +q / /2 ^_ q72;k » +q » /2 



(1 - n F>j8 (k - q/2))n FlJ 8(k' + q'/2)n Fi/3 (k + q/2)£ k _ q/2ik " +q " 7 A +q/2k <_ q < /2 <^' +q ' /2 ,k"-q"/a] ( 36 ) 



In terms of the Bosc fields we have. 



J = A£(q)A£„(-q")r-£? (k',q) 



(37) 



This leads to the answers for T\ and T 2 given in the main text. The arguments in the bose case is analogous. 
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